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The temporal evolution of an atomic wave packet interact- 
ing with object and reference electromagnetic waves is inves- 
tigated beyond the weak perturbation of the initial state. It is 
shown that the diffraction of an ultracold atomic beam by the 
inhomogeneous laser field can be interpreted as if the beam 
passes through a three-dimensional hologram, whose thick- 
ness is proportional to the interaction time. It is found that 
the diffraction efficiency of such a hologram may reach 100% 
and is determined by the duration of laser pulses. On this 
basis a method for reconstruction of the object image with 
matter waves is offered. 
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I. INTRODUCTION 

The success of the last decade in the field of laser 
light cooling below the recoil limit has opened a 

new chapter of atom optics, whose objective is to ma- 
nipulate atomic beams in a way similar to conventional 
optics by exploiting the wave properties of the particles. 
Indeed, if the momenta of cooled atoms approach those 
of photons, diffraction effects may manifest themselves 
particularly strongly during atomic interaction with spa- 
tially inhomogeneous radiation. For the corresponding 
part of the de Broglie wave spectrum, this provides the 
possibility of supplementing the traditional atom-optics 
set of elements such as mirrors ||, diffraction gratings 
J|,[|, or lenses |(| with holograms of different objects, 
the conventional optics analogs of which have been well 
known for several decades j?j . Creation of matter waves 
with the required amplitude and phase characteristics is 
the main task assigned to such atomic holograms. Since 
these characteristics are the same as those of the object 
wave, one obtains a powerful and convenient tool for holo- 
graphic imaging with atoms. The latter may have useful 
practical applications ranging from atom lithography |^] 
to the manufacturing of microstructures or quantum mi- 
crofabrication. 

One possibility to make an atomic hologram is to cre- 
ate a mechanical mask with appropriate transparency 
for the incident atomic beam (the analog of a two- 
dimensional optical hologram). Such a hologram has 
the advantage of being permanent. However, up to now 
only masks with binary transparency have been prepared. 
For example, in the experiment JsJ the mask was written 
onto a thin silicon nitride membrane and allowed for ei- 
ther complete or vanishing transmission of the beam at a 
given point. Since gradually varying transmission of the 



beam is required for the correct holographic storage of 
information, this reduces resolution in the reconstructed 
image. 

A very interesting proposal has recently been reported 
in Ref. ]10|] , where the authors regard a Bose- Einstein 
condensate (BEC) as the registration medium for the 
atomic hologram. In the suggested method the informa- 
tion issuing from an object is encoded into the condensate 
in the form of density modulations by using the object 
and reference laser beams which form the writing optical 
potential. The reconstruction of matter wave arises due 
to s-wave scattering of the reading-beam atoms on con- 
densate inhomogeneities. This proposal illustrates the 
wide potential applicability of the BEC, which, since it 
was realized experimentally pd[ , is now available almost 
routinely in several laboratories. 

In a previous paper p"2] ] we have shown that an atomic 
hologram may also be constructed under certain condi- 
tions as a superposition of reference and object electro- 
magnetic waves, which is common for optical holography. 
The creation of the intended matter wave occurs when an 
ultracold atomic beam is diffracted from this hologram, 
which in fact is just an inhomogeneous light field. The 
main advantages of the proposed scheme are its simplic- 
ity because of bypassing the recording process and, as 
a consequence, the absence of aberrations in the stored 
information. In this sense our approach is close to the 
non holographic scheme of wave front engineering |13|, in 
which the center-of-mass wave function of an atom can 
be arbitrarily shaped by means of a sequence of suitably 
formed laser pulses. 

The main assumption employed in our holographic 
scheme is the linear response of the atomic system to the 
laser-field inhomogeneity. It requires, in particular, that 
the perturbation of the incident atomic beam is weak, 
and sets an upper limit on the object wave amplitude [see 
Eq. (50) in Ref. As a result, only a small portion 

of atoms in the beam can be transferred into the recon- 
structed matter wave. Operating in the linear-response 
regime decreases the diffraction efficiency of an atomic 
hologram, i.e., the ratio of the intensity of diffracted 
atomic waves to the intensity of the reading beam, which 
may be crucial for practical applications. In conventional 
optical holography such a situation corresponds to the 
kinematical regime of information recording . On the 
other hand, the coupled wave theory of Kogelnik Jli) , 
and the theories based on the dynamical approximation 
]l5| p"7| |, provide a recipe for creating a hologram with 
high (up to 100%) diffraction efficiency. To achieve this 
goal it is necessary to control, among other parameters, 
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the weakly inhomogeneous magnetic field onto the path 
of prepolarized particles and appropriately adjusting the 
field gradient, it is possible to suspend the ground-state 
atoms everywhere except in the region of interaction with 
radiation. But if the laser frequency is far from all atomic 
transitions, the contribution to the total force induced by 
spatially dependent shifts of the Zeeman levels is negligi- 
ble. Under this condition, atoms move like free particles 
that are affected only by the electromagnetic waves. 

In Sec. H we specify the magnetic field to compensate 
for the gravity effects and derive an equation that de- 
scribes the dynamics of ground-state atoms interacting 
with the object and reference beams. An approximate 
solution of this equation is found without assuming lin- 
ear response to laser-field inhomogeneity, and its domain 
of validity is determined. For reasonable experimental 
conditions the solution admits an ato m-optics interpre- 
tation, which is presented in Sec. Ill, namely, the inho- 



FIG. 1. Typical layout design of laser beams and matter 
wave packets suitable for atomic holography. 



mogeneous laser radiation is shown to behave like a three- 
dimensional hologram with respect to the impinging wave 
packets. A numerical simulation of such a hologram cre- 
ated with a 31-mode object beam is carried out, and high 
diffraction efficiency is explicitly demonstrated. Section 
[V concludes with a summary of the results. Certain 
mathematical details concerning derivation of the basic 
formulas are relegated to the Appendixes. 



II. BASIC FORMULAS 



A. Compensation for gravity 



the thickness of the hologram. Unfortunately, thick- 
ness controlling is difficult in a scheme of atom hologra- 
phy without a registration medium like the one in Ref. 
@ (see Fig. |l| for details). 

The purpose of the present paper is to suggest an ap- 
proach for creation of atom optics holograms that will 
combine the advantages of our previous scheme with high 
diffraction efficiencies. We will show that high diffraction 
efficiencies can be realized if one restricts the extent of the 
atomic hologram in the time domain rather than in space. 
The hologram will thus work in a pulsed regime pumping 
atoms from the beam or from the initial wave packet into 
the reconstructed wave. Note that the suggested regime 
is well compatible with Raman cooling methods p] (in- 



cluding laser cooling below the gravitational limit |18|) 
and the recent realization of an atom laser Jl9| , which 
are capable of repeatedly reproducing the coherent or al- 
most coherent atomic wave packets necessary for actual 
implementation of a reading beam. 

Another important prerequisite for successful wave 
front reconstruction with massive particles concerns the 
need to compensate for the potentially detrimental in- 
fluence of gravitational effects. Fortunately, the bulk 
of the atoms possess a magnetic moment, and all one 
has to do is use the Stern-Gerlach effect. Superimposing 



To be specific, let us consider an atom with a J = | 
to J = | transition, e.g., sodium or cesium. A magnetic 
field B(r) is applied to compensate for gravity. It is sup- 
posed to contain a homogeneous component Bo directed 
along the gravity acceleration Bo TT g- The remain- 
ing inhomogeneous part of the field Bi(r) = B(r) — Bo 
should be small compared to this component, 



|B 1 (r)|«i?o = |Bo 



(2.1) 



As we will see below, to fulfil this condition it is neces- 
sary to take Bo in the range 10 3 — 10 4 G. In practice 
such a field is strong enough to induce Zeeman shifts 
which considerably exceed the hyperfine splitting inter- 
vals ~ ?iu;hfs (but not the multiplet ones). There- 
fore an internal atomic eigenstate \ J,I,Mj,mi) may be 
parametrized by the set of quantum numbers consisting 
of the angular momenta of the electronic shell J and nu- 
cleus /, and their local projections onto the direction of 
the magnetic field, Mj and mj. The corresponding en- 
ergy eigenvalue is determined not only by the multiplet 
level Ej but also by the magnetic field B(r) = |B(r)| and 
therefore is spatially dependent, 

E\j,i,Mj,mi)( r ) = E J + aMjmj 

+ ( f i B g L Mj- fi nncmi )B(r), (2.2) 
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where a is the hyperfine coupling constant (a oc ?Lu;hfs, 
e.g., for Na a/h = 885.8 MHz), gi denotes the Lande 
factor, and /i nuc is the nuc lear magnetic moment. Be- 
cause of the condition ([O]), such a spatial dependence, 
however, mainly arises from the longitudinal (r) = 
Bo • Bi(r)/_Bo], rather than the transverse [B^r)] com- 
ponent of the vector Bi(r), provided that the compo- 
nents are defined relative to Bo- This is evident from the 
expression 



E a (r,t) 



B(r) 



B, + B\{v) +[Bf( 



B + M'(r) + [B^(r)] /(2B ), 



(2.3) 



where the term containing B^ (r) is small and can be 
neglected. Consequently, by adjusting the gradient of 
the field b\ (r) one can achieve translation invariance of 
the ground state \g) = |l/2,/, —1/2,1) (or another state 
with J = 1/2) in three dimensions: 



E\ g j (r) — Mg • r = const. 



(2.4) 



For example, to balance the gravitational force in this 
way for sodium it is necessary to create a gradient 
V5{(r) = &ig/|g|, where b x = -4.033 G/cm. This 
condition does not contradict the Maxwell equation V • 
Bi(r) = 0, because variation of B^(r) is not restricted. 
Note also that the choice Bq = 10 3 — 10 4 G maintains 
condition (2T) very well within a spatial region of size 
~ 10 cm. 

All the other levels are affected by the residual external 
potential. In particular, the force f e acting on the atoms 
in the excited state, e.g., |e) = |3/ 2, I, —3/2, I), may be 
estimated from Eqs. (2.2) and (2.4) as |f e | ~ Mg. 



B. Dynamics of the ground-state atoms 

In our scheme, we use pulses of laser light at frequency 
ui which is roughly tuned to the \g) — > |e) transition. 
If the typical size 2L of the atomic sample is restricted 
by the condition L <C a/(Mg), one may regard Ei e \(r) 
as the excited level that is closest to resonance within 
the whole interaction domain. Indeed, the maximal spa- 
tial shift of the level ~ MgL induced by the force f e 
appears to be much less than the hyperfine splitting in- 
tervals (MgL <C a ~ 7wa>hfs)j an d the hierarchy of de- 
tunings is retained. Therefore an atom initially in the \g) 
state behaves as a two-level system with respect to the 
processes with stimulated emission of photons. 

The atom moves inside a superposition of the reference 
and the object beams during the whole laser pulse. Each 
beam is represented as a discrete sum of plane monochro- 
matic electromagnetic waves. In particular, we use the 
following decomposition of the electric field in the object 
beam: 



m>l 



E m exp(ik 

m 

r - iwt) + c.c. 



(2.5) 



where E m and k m stand for the complex amplitude of 
the mode to and its wave vector, respectively. Such an 
approach does not restrict the gen erality of our consider- 
ation, because the expression (2J3) must well describe the 
real laser field only in the atom-laser interaction region. 
Evidently, the latter requirement can always be satisfied 
by decreasing the minimal angle between the mode wave 
vectors. In this case we can also regard the reference 
beam as a single mode (with the index m = 0), 



E r (r, t) = E exp(ik • r — iut) + c.c. 



(2.6) 



which is a typical arrangement for optical holography. 

Since the atomic dipole momentum operator d is di- 
agonal in quantum numbers / and mi, the transitions 
that change mj are allowed only due to hyperfine inter- 
action. As a consequence, the excited state |e) decays to 
the lower ones preferentially in the channel |e) — » \g) 
(with the rate 7). This circumstance makes it possi- 
ble to deal with an atom as a two-level system even if 
spontaneous photon emission takes place. However, to 
simplify the consideration, the coherent scattering pro- 
cesses are assumed to dominate the spontaneous emis- 
sion, i.e., the regime |A| ^ 7 is maintained [^0|,^l), where 
A = ui + [_E| 9 }(0) — E^ e ^(0)]/h is the detuning from res- 
onance in the center of the atom-laser interaction region 
(r = 0). Under such a condition the one-particle den- 
sity matrix in the momentum representation [ p2| has an 
obvious time evolution, 

Pab(PuP2,t) = J dp[ J dp2^Gaa'(Pl,Pi,i) 

a'b' 

x G* w (p 2 , p' 2 , t)p a , v (p'i , P 2 , t = 0). (2.7) 

Here indices a,b . . . span the internal atomic states (e, g), 
and G aa '(pi, p' 1; t) is the Green function of the two- 
component Schrodinger equation describing ato mic dy- 
namics during the \g) <-» |e) transitions [see Eq. (Al) for 
the details]. 

For the situation at hand, the upper electronic state 
can be adiabatically eliminated from consideration (see 
Appendix |A[) provided that the detuning A is large 
enough J|,|(J, 



|A| » \n m \,\f e \L/h, 



(2.8) 



where ft m = (e|d-E m |g)/7i is the Rabi frequency of mode 
m. As a result, the dynamics of the ground atomic state 
is completely determined by the equation for the center- 
of-mass wave function ^ s (p, t), 



9 , , s 



Hp) 

+ E 



m>l v n>l 



7 I 



K),t] 
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+9mfpg[P ~ Hk m ~ k ), t] 

+9m^g[P + M k ™ - k ), t] ) , 



(2.9) 



where w(p) denotes the kinetic energy (in units of h), and 
foi fmn, and g m stand for the effective Rabi frequencies, 
introduced by Eqs. (A3). 



C. Evolution of wave packets 

It is known from the theory of thick optical holograms 
that reconstruction of the original (conjugate) object 
wave arises only if the reading beam is directed along (op- 
posite) the reference wave and has the same wavelength. 
In analogy with conventional optics, let us consider for 
definiteness the evolution of an atomic wave packet whose 
spectrum is initially concentrated around the mean mo- 
mentum of photons in the reference beam. In such a 
case one can expect creation of a matter wave similar to 
the forward object wave. Therefore it is convenient to 



seek the solution of Eq. (2.9) as a sum of wave packets 
approaching the plane modes of the hologram I&6 , 

1p g (p,t) = Vm(p-WCn.,*)- (2-10) 



Initially there are no wave packets corresponding to 
the object beam, so that 

lMp,* = 0) = 0, m>l, (2.11) 

and as a consequence, 

^ (P - hk ,t = 0) = 1> g (p, t = 0) = i) g {p). (2.12) 

Population of these atomic motional states (m > 1) arises 
due to coupling with ^o(P)*); the wave packet corre- 
sponding to the reference beam, 



d 



(2.13) 



where 



w m (p) = w(p + 7ik„ 



Depletion of the state with m 
equation 



A + / Q . (2.14) 
= is governed by the 



d 

igl^oiP, t) = W (p)lp (P, tf+Y 9m^m(P. *) + X(P, t), 

m>l 



(2.15) 



which one can obta in after substituting Eqs. (2.10) and 
Q2.13D into Eq. (&M 



Equ ation ( 2.9|) is thus reduced to the system of equa- 
tions (2.13),(2.15[). The advantage of this step becomes 



obvious after making the following self-consistent as- 
sumption about the momentum spectrum of tp m (p,t), 
to > 0, the validit y o f which was verified for the two- 
mode case in Ref. |fL8| . We will suppose below that all 
nonvanishing functions have narrow distributions around 
p = and, as a result, do not overlap in the expression 
for x(p,*), 

x(P,t) = Y (X 9nip m [p + ft( 2k o k„ -k m ),t] 

m>l ^n>l 

+ 9*n^n[P + HKn ~ K),t] 

n>0 

+ 22 fmn 

n>l l>0 

x^[p + h(k -k t -k m +k n ),t] ). (2.16) 



Under these conditions different parts of x(p, t) give 
incoherent contributions, which are small at low g m 
and can be taken into account perturbatively. In zero- 
order a ppro ximation one omits x(p, t) so that the system 
( |2.13 ),( 2T5| ) becomes homomorphic with the rate equa- 
tions describing an (to + l)-level atom. Note that the 
stationary solutions of this truncated system exactly co- 
incide with eigenmodes of the corresponding optical holo- 
gram |16[j . 

To make further progress it is convenient to perform 
the Laplace transformation (to > 0), 



ipm (P, A) = / die A Vm(p,t) 
Jo 



(2.17) 



with the initial conditions (2.11) and (2.12). Then the 



equations for the Laplace transforms will allow an easy 
zero-order solution, 



T(P,A) 



^ ff (p + 7iko), 



(2.18a) 



Vi 0) (P,A)= - ,J? m - #(P,A) ; to>1, (2.18b) 



w m (p) - iX 



where 



r(p,A) = ^( P )-«A+x: 9 r -x - 



(2.19) 



Similarly, the next iteration reproduces the first-order 
solution, 



^ i) ( P ,A)=vr(p,A) + 



(0), 



-X (0) (P,A) 
T(p,X) ■ 



(2.20a) 



^ } (P,A) 



^^4 1} (p,A), to>1, (2.20b) 



w m (p) - «A 
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FIG. 2. Definitions of the values A&, 
simplest case of a two-mode object wave. 



and dp in the 



where X^(p, A) is obtained from the expression ( 2.16| ) 
after making the substitutions ip m (p,t) — > ipmHp, A), 
m > 0. In principle, we can get the solution with any 
desired accuracy by repeating the iterations but it will be 
sufficient to restrict ourselves to the first-order formulas 
in what follows. 

The time-dependent wave functions are obtained as the 
inverse Laplace transforms of "0m (p, A) according to the 
Mellin formula 



ipm(p,t) = 2iri J 



e+ioo 



dAe A Vm(p,A), e>0. (2.21) 



Finally using Eqs. ( 2.2C ) and ( 2.21 ) one can easily obtain 
the ground-state component G 00 (p, p' , t) of the Green 
function appearing in the formula fl2.7|) . The correspond- 
ing expression is derived in Appendix 0. 



D. Validity of assumptions 



admit an approximate analytical solution presented in 
Appendix [ 



The main feature of the near-resonance solution (CI) 
is that the initial atomic wave packet transforms into 
motional states with m > 1 at a time r„ (time of the nir 
pulse H) 



2.9s 



(2n-l), neAf. 



(2.22) 



Here gs stands for the overall effective Rabi frequency 
defined in Eq. flC6|). This transition is velocity selective 
and is most efficient when the Bragg resonance condi- 
tion p • Afe = is satisfied, where A& denotes the typical 
difference between wave vectors in the object and the ref- 
erence beams (cf. Ref. 20 1). The width of the peak in the 
momentum distribution along the direction of vector Afe 
(the interval from the maximum to the first minimum) 
depends on the interaction time, and for t < 2ti is 



S p (t) 



gMgs 
A fe 



(2.23) 



For a given value Afe = | Afc| it decreases with g^. There- 
fore the smaller the effective Rabi frequencies g m the nar- 
rower the momentum spectrum of ipm (p, t). 

To prevent all nonvanishing functions comprising the 
term x(Pj^) horn overlapping in momentum space their 
spectra must be concentrated within the domain |p| < 
hdk at t ~ ri, where 



4 = min |k m - k„| 

m,n>0 



(2.24) 



is the minimal distance between different wave vectors 
of the laser beams (see Fig. ||) . Since the spectral extent 
along the direction of vector Afe is characterized by 5 p (t), 
we immediately get the condition 



Spin) < hS k - 



(2.25) 



In agreement with Eq. (2.23) it sets an upper limit on 
the overall effective Rabi frequency, 



9T. 



< M fc Afe/(2V3M). 



(2.26) 



In the transverse direction the spectra are the same as 
that of the initial wave packet ipg(P + ftko)- Therefore 
one must impose another condition, 



(p' - hk ) x Afc| < MfeAj 



(2.27) 



Green function G gs (p,p',t) 



Let us first check that the zero-order solution fl2.1S| ) which restricts allowed values of p' in the domain of the 
indeed has a narrow momentum spectrum around p = 0, 
provided the initial conditions are chosen properly and 
the effective Rabi frequency g m is small enough. We 
may restrict ourselves to examination of the region D = 



When the inequalities (2.25) and (2.27) are satisfied 



the main correction to the zero-order solution ipm* (p, t) 



{p : [w m (p) — w (p)] 2 < \g m \ , Vto}, where all the func- 



tions in the truncated system of equations ( 2.1 3| ) , ( 2 . 1 E ) 
can come into resonance. In this region the equations 



caused by the term x(P> t) arises outside the near- 
resonance region T> and depends on the geometry of laser 
beams. So for t ~ t\ and a two-dimensional (2D) holo- 
graphic setup like that in Fig. || (i.e., all k m are coplanar 
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FIG. 3. Absolute values of zero-order solution lipg ' 1 ] ( a ) 
and first-order correction to it (b) as functions of the 
momentum component p x and the angle 6 between ki and 
k2. The arbitrary units scale is the same in both plots. The 
remaining components of p are fixed as follows: p y = hko and 
p 2 = 0. The geometry of the laser beams is chosen as in Fig. 
^| The effective Rabi frequencies gi = g% = 10 Hz, /12 =0.1 
Hz. 



vectors) the relative correction has the order of magni- 
tude e r = 8 p (ti) I (h8k) <C 1, as follows from Eqs. ( 2.16 ) 
and (Clb). To obtain this estimate one must apply the 



standard holographic restriction on the intensities of laser 
beams | Eo | 2 3> |E m | 2 , m > 1, which leads to the inequal- 
i ty \g m \ S> |/mn|< Then the third term in the expression 
(2.16) should be discarded, because it is proportional to 
f mn and, consequently, is much less than the first and 
second ones (oc |<7 m |). 

To illustrate the consistency of our approach in the 
case of a two-mode object wave let us consider absolute 
values of the zero-order solution \ipg\ and the first-order 

correction to it, 8^ = \ipg — ^4 l> as functions of the 
momentum component p x and the angle 9 between ki 
and k2. The Cartesian coordinate system is introduced 



FIG. 4. Time dependences of 5^ (solid line), \ipg \ (long 
dashed line), and |^>g | (short dashed line) for a two- mode 
object beam with ki,a -L ko. The components of p are fixed 
as follows: p y — hko and p x = p z = 0. The coordinate axes 
and other parameters are the same as in Fig. |[ 

in momentum space, p = (p x ,p y ,p z ), with the x (y) 
axis directed along (opposite) the vector k (ki). Fig- 
ure U shows the corresponding dependences after 7r-pulse 
time calculated for sodium atoms, provided that the 
initial wave packet has the Gaussian profile ip g {p) oc 
exp [— L 2 (jp — po) 2 /(27i 2 )] with mean momentum po = 
Kko, I ko I = ko — 1.07 x 10 5 cm -1 , and spatial extension 
2L = 0.4 cm. The peaks in the central region of each 
plot correspond to forbidden values of 9 oc 5k < 8 p (ti). 
Outside these peaks (\9\ > 10~ 4 ) the relative correction 
goes down, approaching 0.15 at large 9, that is, below its 
estimation value e r ~ 0.5. 

In the worst case, e.g., k m _L ko, Vm > 1, the estima- 
tion for e r is higher than 8 p (ti) / (TiSk), because it depends 
on 8ui = fi5f./{2M), the minimal kinetic energy an atom 
can get due to transition between the laser modes: 



J sin 2 (g s r) exp(iS u T)d,T 



(2.28) 



Nevertheless, the term x(P> ma V still be treated as a 
small perturbation if the overall effective Rabi frequency 
satisfies the condition, more restrictive than Eq. ( 2.26| ), 



gs < 8a. 



(2.29) 



Otherwise, if g^ ^ 8^, the interaction time should be 
limited so that t <C T\ . Figure || shows time dependences 
of | - 0g°' ) |7 an d \ipg | for a two-mode object beam with 
8 P {i~i)/{h8k) =0.1 and gz = 748^. We see that in the 
considered unfavorable configuration e r does not exceed 
0.1 even if t = 0.5ti. 

In most practical cases, however, only a small fraction 
of las er modes has a geometr y lea ding to the condition 
( |2.29 ), and the requirement ( 2.26 ) appears to be suffi- 
cient. 
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III. ATOM-OPTICS INTERPRETATION 



CT = yjL 2 +ih(t + T)/M, 



(3.6) 



A. General consideration 

In an idealized situation, one may imagine that all 
atoms are initially in a pure state described by the Gaus- 
sian profile 



s(p;po) 



L 3/2 



£3/2^/4 



exp 



-£ 2 (p-Po) 

2h 2 



h p - ro 



(3.1) 



with mean momentum p close to ftk , space position 
r , and ver y sm all dispersion [L 3> %/8 p {t\)\. Accord- 
ing to Eq. (2/7), after interaction with the laser beams 
over a time period r < t\ and subsequent free propaga- 
tion during time t the atoms remain in a pure state, and 
their wave function can be represented as a superposi- 
tion of useful signals ip^'^ (p, r, i; po) and a background 
"0 (b) (P> T >*;Po), where 

^ (CT) (p,r,t; Po ) = e -™(p)*y d P 'G^( P ,p',r). 9 (p'; Po ), 

(3.2) 

a e {s,r,b}. The functions G^ CT )( P ,p',r) are defined by 
Eqs. (B2). In the case considered here they have explicit 



analytical expressions relying on the analogy with Eqs. 
(CI). These expressions are exact in the limit L — > oo, 
and k — * 0, where k — po/h — k> Omitting the ir- 
relevant common phase factor exp[— iko • i"o — ifo T ~ 
{i/K)E\ g \ (0)r] we can readily infer that the Fourier trans- 
form of V^HPi r i t\ Po); 

Vv(r, r, t; po) = A r {%K, r) 7o (r, K y^-r-i W {m )t ^ (3 _ 3) 



propagates like the reference beam. Otherwise, the trans- 
form of ^ (s) (p, r,i;po), 



^ s (r, r,t;p ) = A s {%k,t) 



x^ 7m (r,K)^e ik — <*<•»>*, (3.4) 

generates a matter wave, which inherits the amplitude 
and phase characteristics of the object beam because 
g m oc E m as follows from Eq. (A3c) and the definition of 
the Rabi frequencies f2 m . The last assertion also takes 
into account that all functions 



7m( r >«) 



7r 3 /% 3 



exp 



2 a 



(3.5) 



used in Eqs. ( ft .3] ) and (3.4) slowly depend on r within 
spatial regions ~ 2\a\ 2 / L, each centered around the point 
r™ = r™ + h(t + t)k/M, where 



r " = r 



h(k + k) hk 

Tr-rr T H -— t, 



2M 



71/ 



(3.7) 



and k is some typical wave vector in the object beam. 
The small phase shifts <5^(r — r™) 



Mr) 



2cr 



1 f L * r . K+ Ht + r)(r>-LW )] 



M 



introduced by 7m(r, *c) vanish at small k and large L. 

If the overall effective Rabi fr eque ncy is chosen in 
agreement with the results of Sec. [ID , the background, 
which itself represents a first-order correction to the wave 
function ipg ^ (p, r), appears to be small at any time 



t < T\. Consequently, since the states (3.3) and (3.4) 



are spatially separated after the free propagation period 
tmin — 2LM/(hAk), one may observe a matter wave 
ip s (r, t, t; po) cloning the object beam in the space-time 
region S = {(r,t) : |r - f™| < L, Vm > 1; t > t min }, 
where all atomic wave packets related to different modes 
of this beam still overlap each other. It should be noted 
that S ^ only when the observation time is limited 
by the value i max = LM/[Kksm(6 max /2)], where 6> max 
characterizes the maximal divergence angle of the object 
beam, and k = |k|. In a given context, the physical 
meaning of conditions ( 2.26| ) and (2.29) consists in the 
requirement that a more delicate mechanism (lower laser 
intensity) has to be used in order to restore more de- 
tailed information. The conditions have a counterpart in 
the theory of optical holograms [see, e.g., Eq. (4) in Ref. 
Jl6| ] which, in turn, is responsible for the low intensity 
of noise in the reconstructed wave. 

In a more realistic case we may expect the initial 
atomic state to be a statistical mixture described by the 
density matrix 

P ffs (pi,p 2) 0) = yV/(p')s(pi;pV(P2;p'), (3.9) 

where /(p) denotes a momentum distribution fun ction . 
If this function is compatible with the condition ( |2.27 ), 
one can readily obtain an expression for p gg (pi, p 2 , t) at 
any time. In the region S, it takes the following form in 
the coordinate representation 

pgg(ri,r 2 ,T,t) = J dp7(p')^(ri,r,t;p')C( r 2,T,t;p'). 

(3.10) 



Since A s (p,r) is a sharply peaked function having a 
width S p (t) alo ng th e vector [see Eq. (2.23)], the 
integral in Eq. (3.10) is limited in this direction. Let 
us assume that integration in the transverse directions 
is also restricted to a small domain ~ S p (t) due to th e 
finite spectral width of /(p). Then analyzing Eq. ( 3.1C ) 
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FIG. 5. Diffraction efficiency r\ of atomic hologram as a 
function of time domain r (in units of n pulse) and dimen- 
sionless parameter f. 

in the region |ri — r2 1 <C H/6 p (t) one finds that un- 
der the condition t <C t co h = M/[ksm(6 max /2)S p (T)] the 
density matrix factorizes as a product of coherent states, 



(r,r,t) = C 1 / 2 (r) 



E 



7 m (r,0)^e 



ff™ ik m -r-iu)(fik m )t 

(3.11) 



where 



C(t) = J dp \A s (p, t)\ 2 f(p + Kko). (3.12) 

When deriving the formula ( |3.11 ) we allow for small val- 
ues of phase differences |<^(ri — r™) — (5^(r 2 — r")| <C 
7r, Vm,n > 1, appearing in the integration region at 
t *C t C oh, whence 7 m (r, k) ~ 7 m (r, 0). Note that mutual 
compatibility of time conditions t m - ln < t <C t co h restricts 
the possible structure of the object beam: 



sin(0 max /2) < 



2LkS p (T) 



(3.13) 



We see that the superposition of laser beams selectively 
acts only on those wave packets in the initial representa- 



tion of the atomic density matrix Eq. (3.9) whose spectra 



are concentrated ne ar th e vector ?iko, and thus restores 
the pure state Eq. ( 3.11 ). Therefore the inhomogeneous 
laser radiation proves to behave like a three-dimensional 
hologram with respect to the incident atomic beam (im- 
pinging wave packets). 

One can further establish a close relation between 
an atomic hologram created in a time domain r and 
a permanent optical hologram with the thickness d T = 
Hkor/M along the direction of the reading beam. In- 
deed, as is known from optics, the passage of the reading 
beam through a three-dimensional hologram can be in- 
terpreted as a multiple diffraction in which small waves, 



diffracted from different registration-medium layers with 
equivalent transmission of light, interfere constructively 
to form a reconstructed wave of high intensity. The 
same approach can be used to describe an atom-optics 
hologram as a light structure, inducing an optical poten- 
tial through the atom-laser dipole interaction jl2). Here 
the role of equivalent-transmission layers in the medium 
is performed by the equipotential surfaces. Since d T is 
just the distance the impinging wave packet covers dur- 
ing time t, the numbers of crossed interfaces (layers or 
surfaces) are equal for the atomic and conventional holo- 
grams. Therefore, if there were no difference in the initial 
and boundary conditions, the processes of wave front re- 
construction would be identical in both cases. 

The relation between atomic and optical holograms 
makes it possible to classify atomic holograms as thin 
or thick diffractive optical elements, and use the Talbot 
length LTaibot, i-e., the typical interval between consecu- 
tive interfaces, as a characteristic scale to distinguish be- 
tween the two classes 17 . That is, the hologram can be 



considered as thick (three dimensional) if d T > ^Talbot, 
or in terms of time 



> L Ta \ hot M/(hko). 



(3.14) 



For most holographic setups (for instance, like that in 
Fig. ||) ^Talbot ~ 2n/k . Therefore the criterion (3.14) 
leads to a time domain r larger than the period of atomic 
oscillations. Obviously the latter requirement is well sat- 
isfied for t ~ Ti, the time of the n pulse, provided gs is 
chosen in agreement with the condition (2.26). 



B. Diffraction efficiency 

In a regime where the background is small, we can de- 
fine the diffraction efficiency r] of a hologram as the over- 
all intensity of the modes composing the reconstructed 
wave, provided the initial wave packet is normalized to 
1, 



V(t,Po)= I dp -0 (s) (p,r,i;p o ) 



(3.15) 



It is clear, however, that rj depends on the shape of the 
initial distribution as well. Therefore, to be more spe- 
cific let us assume a Gaussian profile ( |3.l[ ) of the imping- 
ing wave packet with infinitely small dispersion L — * oo. 
Then integration over p' in Eq. (3.2) becomes trivial, so 
that 



77(r,p ) = dp G (s) (p,p ,r) 



(3.16) 



Using the approximate expressions (Clb) and omitting 
negligible interference terms one obtains from the above 
equation 



v(t,Po) = v(t,0 



e + 1 



sin' 



1 



(3.17) 
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where the dimensionless parameter 



(po - fik ) ■ A fc 
M<? E 



(3.18) 



characterizes the deviation of the initial atomic momen- 
tum from the mean momentum of photons in the refer- 
ence beam. 

According to this simple formula the diffraction effi- 
ciency achieves a maximum at r = r n / \J ^ 2 + 1 and can 
reach 100% if f = (see Fig. g). 



C. Numerical example 

In the following we show two-dimensional results ob- 
tained for Na assuming an experimental setup like that 
in Fig. H (i.e., all k m are coplanar vectors and k _L k ). 
The image to be reconstructed is a thin line of width 
A = 27r/fco perpendicular to the laser-beam plane. To 
decrease the amount of computational work we reduced 
the number of object wave modes to 31 and set up 
$max = 7r/4. Such a field approaches the desired sin- 
gle line within a region of size ~ 60A centered aroun d the 
point r = if all the laser modes in the expression (2.5) 



have identical amplitudes E TO , and their wave vectors k n 
are equidistant 



km — kn 



o sin 



7r(m — 16) 



120 



cos 



7r(m — 16) 



120 



, 



(3.19) 



The corresponding profile of the object beam intensity 
distribution I(x) is shown in Fig. where the x axis is 
directed along the vector ko = (ko, 0, 0). The interference 
fringes, which are a consequence of the finite number of 
modes, can easily be separated from the central line and 
therefore do not contaminate our consideration. 

In numerical simulations the optical pulse duration was 
taken to be T\ = 2.82 x 10 -2 s, to demonstrate the high- 
est diffraction efficiency. The remaining laser light pa- 
rameters were fixed as follows: Rabi frequencies VLn = 1 
MHz and TO = 0.01fJ for all 1 < m < 31, detun- 
ing A = -1 GHz (7/ A sa 0.06), the effective Rabi fre- 
quencies g rn = 10 Hz, f mn — 0.1 Hz, and = 55.7 
Hz. Note that for the considered laser-beam geometry 
A fe = V2k = 1.51 x 10 5 cm" 1 and S k = 3.28 x 10 3 
cm -1 , so that the background introduces a relative cor- 
rection of the order e r = 5 p {t\) / (TiSk) = 1.4 x 10~ 2 and 
can be neglected. 

The reconstruction of a real image of the object was 
achieved by impinging Gaussian wave packets ( |3.l| ) hav- 
ing spatial extension 2L = 0.4 cm on the superposition 
of laser beams near the point 



1 
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FIG. 6. Intensity I(x) of the 31-mode object wave as a 
function of observation point r = (a;, 0,0). The inset shows 
optical image of a single line ~ A created in the central region 
~ 60A. 



L 



hkoTi 
M ' tan(# max /2) 



M ' 



(3.20) 



After the interaction with laser radiation is over, these 
wave packets appear at a distance Lj tan(6 l max /2) = 0.48 
cm from the image. As a result, the most intensive mat- 
ter field in the imaging region may be observed after free 
propagation time t = i max cos(6* max /2) = 0.16 s, which 
lies within the limits i m i n = 9.6xl0 _2 s and i max = 0.17 
s. Figure [7] shows the corresponding atomic density pro- 
file p gg (r,r,Ti,t) when the mean momentum of the ini- 
tial wave packet is exactly equal to fi.ko As is seen from 
the bottom part of the plot, the atomic profile displays 
a good match with the distribution of the object beam 
intensity. The at taine d diffraction efficiency calculated 
according to Eq. (3.15) is 98% in this case. 



When initial state is a statistical mixture (3J)) with 
momentum distribution function /(p) uniform along the 
x axis, the atomic density profile acquires a shape rep- 
resented in the Fig. |[ Since condition (3.12) does not 
hold at the chosen laser light parameters, the size of the 
reconstructed line appears to be ~ 4 times wider than 
one might expect from a coherent reading beam. Never- 
theless, such image broadening is not too large, so that 
the atomic hologram can be used even in this unfavorable 
design. 



IV. CONCLUSIONS 

In this paper we have studied a method of driving ul- 
tracold atom propagation using effective holograms made 
of laser radiation in a specified time domain. We have 
shown that the scattered atomic wave packet may inherit 
the features of the object electromagnetic wave provided 
the atomic internal ground state possesses a translation 
invariance due to compensation of gravity with the Stern- 
Gerlach effect. We have established a close relation 








FIG. 7. Atomic density p gg (r, r, Ti, t) as a function of ob- 
servation point r = (x,y,0). The bottom part of the plot 
compares atomic profile (solid lines) with the object beam in- 
tensity distributions I(x) and I(y) (dashed lines) in the planes 
y — and x = 0. 

between the atomic hologram created in the time do- 
main and a thick optical hologram prepared in the cor- 
responding spatial region, and have found a recipe for 
controlling the diffraction efficiency of such an atomic 
hologram by means of varying the time domain. Besides 
adjusting the atom-laser interaction time, another way to 
enhance diffraction efficiency has proved to be the cool- 
ing of the atomic beam so that all the particles get the 
same momentum as the momentum of photons in the ref- 
erence wave. A special role here may be played by BEC 
and coherent atomic-beam generators, which are under 
development now |2S|] . 

We have considered dilute atomic samples, i.e., we have 
not included many-atom interaction s p4| , which may 
lead to nonlinear atom-optics effects p5| along with en- 
hancing the background. The conditions under which 
these interactions can be neglected were elaborated in 
our previous paper ]T^ | using the mean-field approxima- 
tion applied to the Maxwell-Bloch equations |2(j] and are 
well satisfied when the mean-field interaction energy per 
particle is much less than the typical kinetic energy of 
an atom. We have also neglected such possible sources of 
background as spontaneous emission of photons and fluc- 
tuations of the laser frequency. While the first of these 
sources may be eliminated by keeping the laser 



FIG. 8. Atomic density p gg (r,v,T\,t) obtained when the 
initial state is a statistical mixture with uniform momentum 
distribution along the x axis. Other notations are the same 
as in Fig. (?]. 

detuning much bigger than the spontaneous emission 
rate, the second one is determined by the spectral width 
of two-time electromagnetic-field correlation functions 
p2| , p7[ , and substantially decreases if all field modes orig- 
inate from one initial laser mode. 

Although our scheme of an atomic hologram has been 
developed for codirected reading and reference beams it 
can readily be modified for an experimental setup with 
opposite propagation of the beams. In full analogy with 
conventional optics, such a hologram will reconstruct the 
conjugate object wave. 

In conclusion, atom-optics holograms appear to be a 
useful implement for solving some of the basic techno- 
logical problems in the field of atom lithography. For 
instance, they will make it possible to grow 3D circuitry 
components by depositing an arbitrary multilayer pat- 
tern of impurity atoms on a silicon substrate. 

APPENDIX A: ELIMINATION OF THE 
EXCITED STATE 

In the rotating wave approximation the two- 
component Schrodinger equation, rewritten for ground- 
and excited-level wave functions tp g (p,t) and -0 e (p,t) 
slowly varying in time, takes the form 
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[tfl(p)+A]^(p,t) 

-J2 n * m MP + ^ m ,t), (Ala) 

m>0 

Kp)-*fe-V]^ e ( P ,t) 

- J] ^m-0g(P - ^k m ,t), (Alb) 

m>0 



where J7 m stands for the Rabi frequency of mode m 
and the terms u>(p) = p 2 /(2Alh) and — if e • V arise in 
momentum space from the kinetic and potential energy 
(— f e • r) respectively. 

The route by which one can adiabatically eliminate 
the component ip e (p,t) from Eqs. (Al) implies a self- 
consistent assumption \ip e \ <C \ip g \ leading to the zero- 
order solution of Eq. ( |Ala| ): tp g (p,t) ~ cxp{— i[w(p) + 

A]t}ip g (p, t 0). After substitution of this expression 

into Eq. ( Alhj ) the latter is solved perturbatively with 
respect to the potential energy term: 



Vv(p,*) 



\p £l m i>g(p - hk m ,t) 
w(p) ~ w(p - hk m ) ~ A 



+ . . 



(A2) 



where the ellipsis denotes omitted terms that include a 
small (oc |f e |L/|7iA|) first-order correction to ipe(p, t) and 
also terms that oscillate with the nonresonant frequency 
w(p) and therefore give a negligible contribution when 
one u ses the above expression within the context of Eq. 

S- 

For an ultracold atomic sample one can further discard 
the kinetic energy terms in the denominators of the ex- 
pression (|A2|). As a result, Eq. (Ala) takes the form of 
Eq. (2.S), provided that the effective Rabi frequencies /o, 
f mn , and g m are defined as follows: 



Jo-Ty! |^m| 2 , 



A ' 



o o* 



(A3a) 

(A3b) 
(A3c) 



APPENDIX B: GROUND-STATE GREEN 
FUNCTION 

Here we present the first-order approximation to the 
ground-state component of the Green function determin- 
ing the time evoluti on of the atomic density matrix ac- 
cording to formula (2.7): 



G gg (p,p',t) = e(t) G^(p,p',t), (Bl) 

where common phase multiplier e(t) = exp[iwi — 
(i/h)Ei e \(0)t] recovers the solution (2.21) from its slow 
time dependence, and 



G^(p,p',t)=M $\p,p',\) 



G (S) (P,P',*)=E^[^ ) (P.P'^) 

rti > 1 



G (i) (p,P , ,*) = X)M[^(P,P'.A) 

m>0 



(B2a) 



(B2b) 



(B2c) 



In these formulas the operator M. stands for the in- 
verse Laplace transformation and shift of the momentum 
arguments 



M 



( ff )Cr. rV 



P,p', A) 



2m / d\e xt 

J € — ZOO 



x^)(p-^k m ,p',A), (B3) 
e > 0, a £ {0, b}, m > 0, whereas 



CVp',a) 



T(p,X) 



5 3 (p + hk -p'), (B4a) 



T(p,X) 



rX(p,p',A), 



(B4b) 



!ff) (p,p',A) 



,f .^ W,A), (B4c) 
tOm(p) ~ *A 



m > 1, and the expression for x(p, p',A) is obtained 
from Eq. 

X (p, p', A) = ( E [ p + H ( km ~ k »)' P'- A l 

m>l ^ n>0 

+ E ^n^m [P + ^( 2k - k « - km), P', A] 

n>l 

+ E E^ mn 

n>l i>0 

x 0[ h) [p + 7i(ko - k; - k m + k„), p', A] 

(B5) 
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APPENDIX C: NEAR-RESONANCE 
APPROXIMATION 

Here we develop the near-resonance approximation in 
order to get an explicit solution of the truncated system 
of equations (2.13) and (2.15). In the region V one can 



approximately treat the kinetic energy terms w m (p) cor- 
responding to the different modes of the object beam 
(m > 1) as equal: w m (p) w w n (p) sw w(p), where 
w(p) = w(p + Kk) + A + fo, and k is some typical wave 
vector in the object beam. Under this condition the in- 
tegral in Eq. (2.21) can be calculated explicitly, and the 
wave functions tpm^ (p, t) acquire a simple analytical rep- 
resentation, 

4°\p,t) ~ A r (p,t) e -' ib (P)*V> 9 (p + ?ik ), (Cla) 



Vi 0) (P,*) c A s (p,t) 9 -^e- ib ^ g (p + Tiko), m > 1. 
3s 



In these formulas, 



A r (p,i) 



ia{p) 
d(p) 



(Clb) 



sin[d(p)t] + cos[d(p)t], (C2a) 



where 



A s (p,t) = -^sm[d(p)t], 
a(p) 

o(p) = [w(p) - w (p)]/2, 

Hp) = «(p) + w (p), 



d(p) = V«(P) 2 + 9l 



and 



1/2 



. m>l 



(C2b) 

(C3) 
(C4) 

(C5) 
(C6) 



stands for the overall effective Rabi frequency. 
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